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Superconductivity in Naa;Co02 • ^/H20 by charge fluctuation 
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A new mechanism for superconductivity in the newly discovered Co-based oxide is proposed 
by using charge fluctuation. A single-band extended Hubbard model on the triangular lattice is 
studied within the random phase approximation, /-wave triplet superconductivity is stabilized 
in the vicinity of charge-density-wave instability, which is in sharp contrast with the square- 
lattice case. The physical origin of the realization of the /-wave triplet state as well as the 
relevance to experiments are discussed. 
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Recently discovered superconductivity in Naa;Co02 
•yH20^'^^ has attracted great interest since it is real- 
ized in a frustrating system. Compared with high-Tc 
cuprates and Ru-based superconductivity, Sr2Ru04, Co 
oxides have various unique and interesting properties. 
Firstly this is the first Co-based oxide superconductivity. 
Secondly Co atoms form two-dimensional triangular lat- 
tice, which is in sharp contrast with the high-Tc cuprates 
or Sr2Ru04. Thirdly there are three bands close to the 
Fermi surface. Recent intensive experiments have shown 
that the superconductivity is unconventional'^"^^ and the 
material has strong correlation.^'^) 

Superconductivity in the systems with frustration is 
a very interesting subject. Originally the resonating 
valence bond (RVB) idea by Anderson was proposed 
in the frustrating triangular Heisenberg spin system.^) 
This concept was developed to the idea of RVB su- 
perconductivity which is expected to be realized in the 
doped Heisenberg system like high-Tc superconductiv- 
ity. Since the high-Tc cuprates are based on the square 
lattice without frustration, superconducting materials in 
frustrating systems have been long desired. The new su- 
perconductor, Naa;Co02-2/H20, is in this sense an ideal 
material which has a perfect triangular lattice of Co 
atoms. Basically frustration destroys the long-range or- 
der like antiferromagnetism and induces large fluctua- 
tions. There is a potential possibility that there appears 
high-Tc superconductivity using these large fluctuations. 
One candidate in the Co oxides is spin fluctuation, which 
, has been discussed using t-J models on the triangular lat- 
■ tice."'^^"^^-' Another candidate is charge fluctuation which 
is realized in the extended Hubbard model. They contain 
new physics induced by frustration. In the present pa- 
per, we focus on the latter possibility in Naa;Co02-j/H20 
system. 

Let us first discuss the effects of multi-orbitals.^^' ^^-^ 
In the case of Sr2Ru04, the three orbitals {dxy, Ayz 
and dzx) do not mix so much. On the contrary, in the 
present material, Na:rCo02, the three orbitals mix rel- 
atively strongly and three bands are formed near the 
Fermi energy. However it has been pointed out that the 
main band centered around the F point [k = (0,0)) is 
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made out of the symmetric combination of the orbitals, 
i.e., a2g-orbital, or -^{dxy + dj,^ -|- d^x)-^^^ Thus, in the 
present paper, we use a single band effective Hamiltonian 
as a first approximation. 

The t-J model on the triangular lattice has been ex- 
tensively studied by various groups. ^^^^^^ However ev- 
ery mean-field-type theory has concluded that the most 
probable RVB superconductivity has dx2-y2-\-\dxy-^sye 
symmetry. This state has been confirmed by high- 
temperature expansion studies. "'^^^ However, recent /iSR 
experiment has shown that no evidence for the bro- 
ken time-reversal-symmetry is found, which contradicts 
with the prediction in the t-J model. Alternatively it has 
been proposed that charge ordering is realized near the 
electron density n = 1-1-1/3 and n = 1 + 2/3 (more-than- 
half-filling).^°) It is apparent that the 1/3- or 2/3-filling 
is a special filling for the triangular lattice. Motrunich 
and Lee'^^) have discussed a strongly correlated electron 
system just next to the charge ordering. 

Having these in mind, we discuss superconductivity 
due to charge fluctuation in the vicinity of charge density 
wave (CDW). We use a single-band extended Hubbard 
model on the triangular lattice which is given by 

H = ^Ekcl^Cka + U ^n^^Tiii +V^ninj , (1) 

ka i (ij) 

where Cj[^ and Cka denote electron creation and anni- 
hilation operators respectively and {ij) represents the 
summation over the nearest neighbor pairs. £k is the 
dispersion relation of the triangular lattice, namely £k = 
~~2t{cos ky -\- 2 cos ^kx cos ^ky) — fi. In this extended 
Hubbard model, we expect a CDW state when V is very 
large and U > V. For the intermediate V, supercon- 
ductivity is predicted in the vicinity of the CDW insta- 
bility. ^^"^^^ Since previous theories only discussed the 
square or cubic lattice, it is necessary to study the case 
of triangular lattice. Surprisingly we find that the triplet 
/-wave state is stabilized, instead of dxy-'wave obtained 
in the square lattice. 

The LDA calculation^^' showed that there is a large 
hole-like Fermi surface around the F point so that the 
effective hopping integral is t < 0. For the Co oxides, 
the electron density n is more than unity {n is around 
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1.3). In the following, wc use a convention with t > and 
n < 1, which is equivalent to the case with t < and 
n > 1 by electron-hole transformation, Ck — > c]^- The 
"hole-doped case" in the present paper corresponds to 
the Co oxides. We also study the "electron-doped case" 
for comparison. 

To study the effective interaction between electrons 
arising from charge and spin fluctuations, we use random 
phase approximation (RPA).^^"^^) In RPA, the effective 
pairing potentials for the singlet and triplet channels are 
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where V{q) = 2V{cos qy + 2 cos cos ^qy) and lui is 
the Matsubara frequency. Xs and Xc are spin and charge 
susceptibilities, respectively. These quantities are calcu- 
lated within RPA as follows. 



Xs(9,w;) = xo(g,w;)/[i - f/xo(g,wO], 

Xc{q,LOi) ^ Xo{q,L^i)/[l + {U + 2V{q)}xoiq,oJi)]. 
Here xo is the bare susceptibility given by 
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Note that the terms proportional to Xc in eqs. (2) and (3) 
represent effective pairing potentials due to charge fluc- 
tuation. This charge fluctuation contributes equally to 
V and V* since it comes from charge degrees of freedom. 

Figure 1 shows the momentum dependences of the spin 
and charge susceptibilities at T = 0.01. In the following 
we use f as a unit of energy and fee = 1- We calculate 
both hole-doped (n < 1) and electron- doped (n > 1) 
cases, since there is no particle-hole symmetry in the tri- 
angular lattice. Let us discuss the hole-doped case with 
n = 0.8 (Fig. 1(a)), which corresponds to the Co oxides. 
When V = (i.e., {U,V) = (3.64,0)), the spin fluctua- 
tion is much larger than the charge fluctuation because 
the system is close to the spin density wave (SDW) insta- 
bility. Xs has a peak at Q = (0, ^tt) =(K-point) due to 
the nesting condition of the system. On the other hand, 
when V ^ 0, i.e., {U,V) = (3.21,1.2), Xc has a peak 
structure around q = Q where V{q) has its maximum 
value. This shows that V induces the charge density 
wave (CDW) instability for large-F. 

For the electron-doped case {n = 1.2), in contrast, the 
peak in Xs for = moves away from Q as shown in Fig. 
1(b). Near the CDW instability i.e., {U, V) = (2.0, 1.06), 
on the other hand, Xc becomes large at q = Q. Appar- 
ently the peaks of Xs and Xc are located at different posi- 
tions from each other in contrast to the hole-doped case 
where both spin and charge fluctuations become large 
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Fig. 1. Spin (xs) and charge (xc) susceptibilities for (a) the hole 
doped (n = 0.8) and (b) the electron doped (n = 1.2) cases at 
T = 0.01. The parameters are chosen at (U, V) = (3.64, 0) and 
(3.21,1.2) in (a) and {U,V) = (2.96,0) and (2.0,1.06) in (b), 
respectively. 

Next, we calculate the pairing potentials, V and V*, 
as shown in Fig. 2. For ^ = 0, the momentum de- 
pendence of is similar to that of x.5 in both hole- 
and electron-doped cases because contributions from Xc 
is negligible. V* has the same property as although 
the sign is opposite and the magnitude is smaller than 
V^. Near the CDW instability (^ ^ 0), on the contrary, 
both and V* show a negative (i.e., attractive) peak 
at q = Q. Note that constant shift in V* or does 
not affect anisotropic superconductivity. In the region 
away from Q, remains positive due to the contribu- 
tion from the spin fluctuation, whereas V* is negative for 
almost all region because both spin and charge fluctua- 
tions lead to the attractive interaction as shown in eq. 
(3). 

To obtain the onset of the superconducting state, 
we solve the linearized Eliashberg's equation within the 
weak-coupling theory: 

AA(fe) = W-*)(fe-fe',0)^^^^if^A(fe') , (7) 

with Cfc = £fc ^ M- The transition temperature is de- 
termined by the condition, A = 1. In the weak-coupling 
theory, u) dependence of the order parameter A(fc) is ne- 
glected. Although this approximation is quantitatively 
insufflcient, it is expected to be valid for investigating the 
pairing symmetry of A(fe) and for grasping the basic idea 
of the superconductivity mediated by charge fluctuation. 

The point symmetry group of the triangular lattice is 
De, which has six irreducible representations as shown in 
Table I. According to these representations, we classify 
A(fc) or eigenfunctions of Ehashberg's equation. A(fc) 
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Fig. 2. Effective interactions for singlet {V") and triplet (V*) pair- 
ing for (a) hole-doped and (b) electron-doped cases, respectively 
at T = 0.01. 



belonging to Ai, and Ei symmetry correspond to 
spin-singlet states, whereas those belonging to B\, B2 
and El correspond to spin-triplet states. We find that 
A(fe) with maximum eigenvalue belongs to E2 for the 
spin-singlet case and B2 for the spin-triplet case, respec- 
tively. This means that the realized superconductivity 
is d-wave singlet or /-wave triplet. Actually there are 
two degenerate states in the rf-wave singlet channel be- 
cause it is two-dimensional representation (£^2)- These 
two states form d^2_y2-|-ida;y state below predicted in 
the t-J model. Interestingly we did not find f)-wave state 
for the triplet c;liannel in contrast to the Sr2Ru04 case 
where Px+'vpy state is realized. 



Table I. Irreducible representations (IR) of Dq 
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The obtained phase diagrams on the {U, V) plane for 
n = 0.8 and n = 1.2 cases are shown in Figs. 3(a) and 
3(b), respectively. The temperature is fixed at T = 0.01. 
The dashed lines which correspond to the divergence of 
Xs or Xc determine the boundary of SDW or CDW state, 
respectively. 

Let us discuss the hole doped case (n = 0.8) first. In 
Fig. 3(a), the eigenvalue with the E2{d) symmetry be- 
comes larger than unity in the right-hand side of the solid 
line. In this case, the effect of V suppresses the d-wave 



pairing, which is understood as follows. For n = 0.8 and 
without y, the effective interaction for the singlet pair- 
ing is large at q = Q as shown in Fig. 2(a), inducing the 
d-wave superconductivity near the SDW. In the presence 
of V, however, the charge fluctuation gives an opposite 
contribution aX q — Q and thus suppresses the pairing 
potential V^, which results in the suppression of rf-wave 
pairing. 
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Fig. 3. Phase diagram on the {U, V) plane at T = 0.01 for (a) 
hole-doped and (b) electron-doped cases. The dashed lines cor- 
respond to the CDW and SDW instabilities. On the solid and the 
dotted lines, the eigenvalues which belong to E2(rf) and B2(/) 
symmetry reach unity, respectively. In (a), the d-wave pairing 
is dominant near SDW and the /-wave pairing is dominant near 
CDW. In (b), the eigenvalue of the /-wave solution reaches unity 
only in the region where that of the d-wave solution is larger than 
unity. 

In the proximity to the CDW boundary, the supercon- 
ductivity with B2{f) symmetry becomes dominant. The 
momentum dependence of the order parameter is shown 
in Fig. 4. The /-wave solution has three peaks with 
the same sign which arc connected by the wave vectors 
(0, Itt), ("^"■i ~|^) and {:^T^, §''■) in the triangular lat- 
tice Brillouin zone. The order parameter with this prop- 
erty is favorable because of the momentum dependence of 
V* which gives large attractive interactions at these wave 
vectors. In fact, for the triplet pairing, both spin and 
charge fluctuations give attractive interactions and work 
cooperatively since their magnitude are large at the same 
wave vector, q = Q, as shown in Fig. 1(a). Furthermore, 
in the real-space picture, we can understand the stabil- 
ity of the /-wave triplet state as follows. The effect of V 
repels electrons from the nearest-neighbor sites. Then it 
is natural that the amplitude of the order parameter in 
real space becomes large at the six next-nearest-neighbor 
sites. After Fourier transformation, this real-space order 
parameter results in the /-wave symmetry in fc-space. 

Next we consider the electron doped case, (n = 1.2) 
(Fig. 3(b)). In this case, the obtained phase diagram 
differs from the hole-doped case. This difference mainly 
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comes from the momentum depcndcincc; of Xs which has 
peaks away from Q in contrast to the hole-doped case. 
As a consequence, for the singlet pairing, both spin and 
charge fluctuations mediate rf-wave superconductivity as 
shown in Fig. 3(b). As for the triplet case, the effect of V 
is necessary for the stability of the /-wave pairing. The 
peaks in V* in Fig. 2(b) for F = case do not give /- 
wave state because the peak positions are away from Q. 
Near the CDW instability, becomes large but is 
also enhanced. As a result, d- and /-wave states compete 
with each other. 




fluctuation in the vicinity of the ferromagnetic region. 

Nearly ferromagnetic spin fluctuation has been observed 
in Co-NQR experiments. The other mechanism can 
be in the weak-coupling theory. Starting from the Hub- 
bard model on the triangular lattice or in a tight-binding 
model, /-wave superconductivity has been discussed us- 
ing the third-order perturbative calculations.^^' In order 
to determine which mechanism is probable, it is neces- 
sary to investigate the charge fluctuation in the material 
experimentally. 

In summary we have shown that /-wave triplet super- 
conductivity is realized in the vicinity of CDW instability 
in the triangular lattice. Superconductivity induced by 
charge fluctuation is a very new and interesting phenom- 
ena. In the square or cubic lattice, rfj;y-wave singlet su- 
perconductivity was discussed so far.^^"^"'' However the 
situation is drastically different in the triangular lattice 
as shown here. 

The reason why the /-wave state is stable is summa- 
rized as follows. 1) Charge fluctuation is equally helps 
singlet and triplet pairing, since it is the charge degrees of 
freedom. Symmetry of superconductivity is determined 
by the geometry of the Fermi surface. 2) Due to the 
effect of the nearest-neighbor repulsion V , the Cooper 
pairs tend to be formed on the next-nearest-neighbor 
sites avoiding the nearest-neighbor sites. 3) Since there 
are six next-nearest-neighbor sites in the triangular lat- 
tice, the sign of the Cooper pairing can take (H 1 1 — ) 

which fits very well to the lattice structure. 

The authors would like to thank T. M. Rice, P. A. Lec. 
B. S. Shastry, G. Baskaran, K. Ishida, K. Yoshimura, Y. 
Suzumura, and Y. Tanaka for very useful discussions. 
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Fig. 4. (a)Thc momentum dependence of the gap function A(k) 
for the /-wave triplet pairing obtained at (U, V) = (3.21, 1.2) 
and T = 0.01. The contour plot is also shown in (b). 

Experimentally it has been proposed that the triplet 
superconductivity is realized in Naa,Co02-yH20 from the 
Knight shift measurements,^' although there is still a 
controversy about the results.^' Furthremore the relax- 
ation rate 1/Ti measurements in Co-NQR^' suggests 
the presence of the line nodes. Prom these experiments 
the /-wave state found in this paper is promising. 

There are other possibilities for the superconduct- 
ing mechanism. One is due to ferromagnetic spin- 
fluctuation. Although the LDA calculation shows that 
the Co oxides has t < and n > 1, Kumar and Shas- 
tryi3) have proposed a model with opposite sign of t 
{t > 0), which corresponds to t < and n < 1 in the 
present notation. Actually the high-temperature expan- 
sion study in the t-J model on the triangular lattice has 
shown that there is a ferromagnetic instability in the 
case oi t < 0.^^' It is possible that the triplet super- 
conductivity is induced by using the ferromagnetic spin- 
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